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LONG  TERM  GOALS 

To  develop  passive  modalities  of  acoustic  ocean  monitoring  techniques  such  as  acoustic  tomography  or 
acoustic  thermometry  using  Cross-correlation  processing  of  ocean  ambient  noise 

OBJECTIVE 

Coherent  processing  of  ambient  noise,  based  on  these  computed  noise  cross-correlation  between 
receivers,  has  been  suggested  as  a  potential  means  for  developing  passive  methods  for  ocean  acoustic 
tomography  or  thermometry  [Roux  et  al.,  2004].  However,  this  technique  can  require  long  recording 
times  to  extract  stable  coherent  arrivals  from  ocean  noise  correlations  [Sabra  et  al.,  2005a;  Godin,  2010]. 
Array  beamforming  can  enhance  the  emergence  of  such  deterministic  arrivals  between  receiver  arrays. 
For  instance,  a  single  vertical  array  was  used  to  improve  passive  imaging  of  seabed  layers  [Siderius  et  al., 
2010;  Siderius,  2011;  Traer  and  Gerstoft,  2011]. 

The  main  objective  of  this  year  research  was  to  develop  a  general  spatiotemporal  filtering  (STF) 
technique  to  enhance  the  emergence  rate  of  coherent  wavefronts  between  two  spatially  separated  receiver 
arrays.  The  performance  of  this  STF  technique  was  demonstrated  using  ocean  ambient  noise  recordings 
dominated  by  non-stationary  shipping  noise  to  monitor  the  arrival-time  fluctuations  of  coherent  noise 
wavefronts  over  a  six  day  long  observation  time. 


WORK  COMPLETED 

Ambient  noise  was  recorded  on  two  vertical  line  arrays  (VLAs)  separated  by  450  m  and  deployed  in 
shallow  water  (depth  -150  m)  off  San  Diego,  CA  continuously  for  six  days.  Recordings  were  dominated 
by  non-stationary  and  non-uniform  broadband  shipping  noise  (250  Hz-1.5  kHz).  Stable  coherent  noise 
wavefronts  were  extracted  from  ambient  noise  correlations  between  the  VLAs  during  all  six  days  by 
mitigating  the  effect  of  discrete  shipping  events  and  using  array  beamforming  with  data-derived  steering 
vectors.  This  procedure  allows  the  tracking  of  arrival-time  variations  of  these  coherent  wavefronts  during 
six  days  and  may  help  developing  future  passive  acoustic  tomography  systems.  A  detailed  summary  of 
the  completed  work  is  provided  hereafter. 
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Fig.  1  (a)  Bathymetry  map  (depth  denoted  in  fathoms)  in  the  vicinity  of  the  harbor  of  San  Diego,  CA.  The 
test  site,  marked  by  a  box,  is  located  on  the  shallow  Coronado  Bank  (water  depth  ~150nt).  The  zoomed-in 
inset  in  the  bottom-right  corner  sketches  the  relative  orientation  of  the  two  VLAs  with  respect  to  the 
shoreline,  (b)  Averaged  spectrogram  of  the  ambient  noise  (across  the  all  16  elements  of  VLA1)  recorded 
on  the  third  day  of  the  experiment  in  the  frequency  band  [250  Hz-12.5  kHz].  The  vertical  axis  represents 
the  acoustic  frequency  in  a  logarithmic  scale.  The  horizontal  axis  indicates  Greenwich  Mean  Time. 

Broadband  ambient  noise  [250  Hz-12.5  kHz]  was  recorded  continuously  for  six  days  from  January  31  to 
February  5,  2009  (denoted  as  Days  1-6)  near  San  Diego,  CA  using  two  identical  VLAs  moored  450  m 
apart  in  a  shallow  and  nearly  range-independent  section  of  the  Coronado  Bank  (water  depth  —150  m,  Fig. 
1(a)).  Each  VLA  had  16  elements  uniformly  spaced  by  1  m;  the  first  element  being  ~7  m  above  the 
seafloor.  Other  technical  features  of  the  hydrophone  array  deployment  and  the  electronic  system  have 
been  described  previously  [Skinner  et  al.,  2008;  Leroy  et  al.,  2012], 

Due  to  the  proximity  of  San  Diego’s  harbor,  acoustic  data  were  dominated  by  non-stationary  shipping 
noise,  with  episodic  loud  discrete  shipping  events,  for  frequencies  below  2  kHz  (see  Fig.  1(b)).  In  order 
to  mitigate  the  effects  of  these  episodic  shipping  events,  the  continuous  noise  recordings  from  each 
hydrophone  were  segmented  in  1  minute  long  intervals  and  then  homogenized  using  the  same  two 
processing  steps  described  in  previous  studies  [Sabra  et  al.,  2005b;  Leroy  et  al.,  2012]:  (1)  whitening  the 
amplitude  spectrum  of  the  data  in  the  most  energetic  frequency  band  [250  Hz  -  1.5  kHz]  to  diminish 
eventual  strong  spectral  peaks  and  (2)  clipping  the  signal  amplitudes  above  a  threshold  equal  to  three 
times  the  average  standard  deviation  of  the  frequency  whitened  time  series. 

3.  Spatial  origin  of  the  coherent  noise  field 


Sj2) 


(t;k) 


S{1)(t •  k) 

Given  each  1  minute  long  frequency-whitened  and  clipped  time-series  '  v  ’  ’  and 
recorded  respectively  by  the  ith  and  jth  hydrophone  of  VLA1  and  VLA2  (i,j=1..16)  during  the  kth  minute  of 
the  six  day  long  recording  (k=1..8640),  their  normalized  cross-correlation  function  C:i(t;k)  was  computed 
using: 


C, (f, k )  =  J S,(1) (r; k)S(p ( T  +  t-k)dz  /  Jj(SP(r,k))2  Jr,Jj(Sj2)(r;£))2 


dr 


(1) 
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Subsequently,  the  ensemble  averaged  cross-correlation  waveform 


RJt-L,N) 


is  defined  hereafter  over 


a  N  minute  long  recording  interval  starting  at  the  L  minute  of  the  six  day  long  recording, 


R,At,L,N)  =  ^Zc,.>m 

Fig.  2(a-d)  displays  the  coherent  wavefronts  extracted  from  the  ensemble  averaged  cross-correlation 

Ri,i  (fi  L,  N)  over  a  22  hour  iong  interval  (i.e.  N=1920)  starting  on  either  the  third  day  (L=2881)  or  fifth 
day  (L=5761)  of  the  experiment  for  both  positive  and  negative  time-delays.  Based  on  Eq.  (1),  the  main 
coherent  wavefront  in  Fig.  2(c-d),  centered  on  the  negative  time-delay  t=  -0.3  s  (consistent  with  the  VLAs 
separation  of  ~450  m)  results  from  coherent  ambient  noise  propagating  successively  between  VLA1  and 
VLA2  (i.e.  along  the  northward  direction,  see  inset  on  Fig.  1(a))  for  Day  3  or  Day  5.  Furthermore,  the 
similarity  of  the  arrival-times  of  these  coherent  wavefronts  obtained  (Fig.  2(c-d))  is  indicative  of  the 
relative  stability  of  the  water  column  between  the  VLAs  (especially  at  depths  below  138  m)  during  this 
experiment  between  Day  3  and  Day  5.  Conversely,  coherent  wavefronts  occurring  at  positive  time  delays 
correspond  to  coherent  ambient  noise  propagating  successively  between  VLA2  and  VLA1  in  the 
southward  direction  (Fig.  2(a-b)).  Due  to  the  proximity  of  the  main  shipping  lanes  exiting  San  Diego’s 
harbor  just  north  of  the  experimental  site  (see  Fig.  1(a)),  it  was  found  that  these  coherent  wavefronts  for 
positive  time-delays  were  more  dominated  by  isolated  loud  shipping  events  while  the  coherent  wavefronts 
for  negative  time-delays  resulted  from  more  diffuse  shipping  noise.  This  explains  the  slightly  different 
wavefront  structures  and  time  delays  of  the  positive  arrival  wavefronts  (Fig.  2(a-b))  when  compared  to  the 
negative  arrival  wavefronts  (Fig.  2(c-d)).  For  instance,  the  coherent  wavefronts  shown  on  Fig.  2(b)  were 
likely  dominated  by  a  single  loud  shipping  event  passing  close  to  the  VLAs  and  thus  generating  a  more 
complex  multipath  structure  when  compared  to  Fig.  2(c-d). 


Fig.  2  Spatial  temporal  representation  (in  logarithmic  scale)  of  the  coherent  wavefronts  obtained  from 


the  ensemble  averaged  cross-correlation  hJ  ’  ’  (see  Eq.  (2))  between  hydrophone  i=4  of  VLA1 
and  all  receivers  of  VLA2  (plotted  by  associated  depth),  for  a  32  hour  (i.e.  N=1920)  long  recording 
interval  starting  at  the  first  minute  of  each  respective  day;  (a)  Day  3  positive  arrivals,  (b)  Day  5  positive 
arrivals,  (c)  Day  3  negative  arrivals,  and  (d)  Day  5  negative  arrivals.  The  maximum  value  in  each  plot 
was  set  to  unity  (i.e.,  0  dB).  (e)  Comparison  of  the  evolution  of  the  peak  SNR  (see  Eq.  (3))  of  the 
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,(t;L,N ) 


ensemble  averaged  correlation  waveform  ‘  ,=4,-/=2  " ’  computed  between  an  arbitrary  selected  pair 
of  hydrophones  of  the  VLAs  for  increasing  averaging  time  (N=l ...  1440)  starting  from  either  the  first 
minute  of  Day  3  or  Day  5.  The  evolution  of  the  peak  SNR  is  computed  for  positive  and  negative  time- 
delays  in  the  same  time  windows  around  the  arrivals  as  shown  on  (a-d).  (f)  Same  as  (e)  using  instead  the 

peak  SNR  of  the  beamformer  output  >  L,  .  The  four  sets  of  coherent  wavefronts  shown  in  (a-d)  were 
to  determine  the  four  different  sets  of  steering  vectors  used  to  compute  the  beamformer  output 
N  )  jor  eacjq  case 


4.  Arrival-times  variations  of  coherent  noise  wavefronts 


A  coherent  arrival  is  clearly  visible  on  the  averaged  cross-correlation  waveform  ,,J  ’  ’ 
only  when  this  arrival’s  amplitude  (the  signal  of  interest,  e.g.  see  coherent  wavefront  on  Fig.  2(a)) 
becomes  larger  than  the  peak  amplitude  of  the  residual  temporal  fluctuations  (which  sets  the  background 


'‘noise”  level  of  this  waveform)  in  J  ^  occurring  around  the  coherent  arrival  time 


-0.3  s. 


Assuming  that  these  residual  temporal  fluctuations  behave  as  Gaussian  random  signals,  their  peak 

amplitude  can  be  estimated  by  three  times  the  standard  deviation  of  |-gencjat  an(j  piersol5 

2000].  Flence,  the  peak  signal-to-noise  ratio  SNR(L,N)  aVeraged  cross-correlation  waveform 

Rij(t;L,N) 


is  defined  as: 

f 


SNR  ( L,N )  =  10  log 


max 

teTw 


Ri  At;  L,N) 


i  i  L+N  (  i 

—  — y  (c,.  ft-k))2  -  \  — 
\  N  y  N  £l  '’j  (  N 


£  Cu(f,k) 


2  A 


'  J 


(3) 


where  the  maximum  amplitude  search  is  performed  in  a  20  ms  long  time  window  Tw  centered  on  the 
expected  arrival  time  |t| — 0.3  s  of  the  main  coherent  wavefronts  (see  Fig.  2).  The  denominator  in  Eq.  (3) 

represents  three  times  the  standard  deviation  of  the  ensemble  averaged  correlation  assuming 

that  the  fluctuations  amongst  the  individual  cross-correlation  functions  Cij(t;k )  are  uncorrelated  and  have 
the  same  standard  deviation  [Weaver  and  Lobkis,  2005]. 

Fig.  2(e)  compares  the  evolution  of  this  peak  value  of  SNR(L,N)  por  positive  and  negative 
time  delays  for  increasing  averaging  duration  up  to  24  hrs  (1<N<1440)  using  noise  data  recorded  on 
either  the  third  (L=2881)  or  fifth  (L=5761)  day  of  the  experiment  between  the  fourth  element  of  VLA1 
and  the  second  element  of  VLA2  (i=4,j=2).  Note  that  none  of  the  peak  SNR  curves  increases 

proportionally  to  the  theoretical  predictions  of  derived  for  a  stationary  and  isotropic  ambient  noise 
field  [Sabra  et  al.,  2005b;  Weaver  and  Lobkis,  2005].  Indeed,  the  peak  SNR  for  both  negative  and 
positive  time-delays  builds  up  non-uniformly  primarily  due  to  the  occurrence  of  loud  shipping  events; 
especially  around  N=3  00-460  min  on  Day  3  for  positive  time  delays  (solid  and  dashed  lines)  or  around 
N=650-750  min  on  Days  3  and  5  for  negative  time  delays  (plus  and  dot-dash  lines).  Additionally,  the 
peak  SNR  remains  low  (<5  dB)  even  with  long  averaging  times  for  negative  time-delays. 

In  order  to  enhance  the  emergence  rate  of  coherent  wavefronts  across  all  six  days  of  the  experiment,  the 
STF  methodology  introduced  by  Leroy  et  al.  (2012)  was  implemented  to  beamform  the  averaged  cross¬ 
correlation  waveforms  (jj=|  ...M=16).  At  the  frequency  f  the  beamformer  output 
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B(J,  L,  N)  C0rrCsp0nding  to  the  Mxl  steering  vectors  ^  and  ^(f)  applied  respectively  to  the 
M=16  elements  of  VLA1  and  VLA2  is  given  by: 


B(f,L,  N)  =  W1h  (/)  R  (f,L,N)  W2  (/) 


(4) 


where  the  superscript  H  denotes  a  complex  transpose  operation  and  is  the  MxM  cross¬ 
covariance  matrix  between  VLA1  and  VLA2.  The  (i,j)th  entry  of  the  matrix  -g  Fourjer 

transform  of  averaged  cross-correlation  waveforms  /  ^  at  the  frequency  f.  In  order  to  assess  the 

performance  of  the  STF  methodology,  the  steering  vectors  ^  ^  ^  and  ^  were  computed  hereafter 

for  four  different  cases.  As  suggested  by  Leroy  et  al.  (2012),  the  steering  vectors  ^  and  for 

Cases  1-4  were  set  as  the  first  left  and  right  singular  vectors  of  four  different  reference  covariance 

matrices  ^  ’  corresponding  to  the  four  sets  of  coherent  wavefronts  shown  on  Fig.  2(a-d).  Each 
reference  is  a  32  hour  long  average  (N=1920),  which  starts  on  Day  3  (L=2881)  or  Day  5  (L=5761)  for 
positive  arrivals  (Case  1  and  Case  2)  or  for  negative  arrivals  (Case  3  and  Case  4). 

The  cross-correlation  of  time  domain  beamforming  on  VLA1  and  VLA2,  ^’^was  obtained  from 

the  inverse  Fourier  transform  of  the  beamformer  L,  N)  compUlct|  Vv  ith  Eq.  (4)  using  each  of  the 

four  sets  of  steering  vectors  ^  and  ^(f)  defined  above  (Cases  1-4).  The  SNR  of  the  peak 

arrival  of  was  defined  in  a  similar  fashion  to  Eq.  (3).  Overall,  compared  to  the  results  for  the 

element-to-element  correlations  shown  in  Fig.  2(e),  Fig.  2(f)  indicates  the  peak  SNR  of  the  beamformer 

B(t,L,N )  grows  more  rapidly  and  more  uniformly  during  the  same  24  hour  averaging  interval.  For 
instance,  the  peak  SNR  of  the  beamformer  exceeded  an  arbitrary  threshold  of  5  dB  using  an  N=480  min 
(i.e.  8  hrs)  recording  window  for  Cases  (1-4). 

Finally,  the  four  different  sets  of  steering  vectors  (Cases  1-4)  were  successively  applied  to  the  entire  six 

day  long  dataset  to  beamform  all  cross-covariance  matrices  -  480)  (i<l<8160)  computed 

using  a  short  moving  average  of  N=480  min.  Fig.  3  displays  the  time  shifts  of  the  peak  coherent  arrival 

of  the  corresponding  beamformer  -  480)  ^por  £ases  j_4)  vs  indices  L  (i.e.  the  beginning  of 

the  moving  time-window  along  the  “Slow  Time”  axis)  when  the  peak  SNR  was  sufficiently  high  (i.e.  here 
above  an  arbitrary  threshold  of  1.5  dB)  to  ensure  clear  arrival  tracking.  Overall,  this  STF  technique  was 
robust  enough  to  extract  a  small  and  quasi-linear  upward  trend  of  the  peak  coherent  arrival  of  the 

beamformer  -  480)  ^ming  the  entire  six  days  of  the  experiment  for  all  four  different  sets  of 

steering  vectors  (see  Fig.  3).  A  comparison  between  the  positive  arrivals  Fig.  3(a-b)  and  the  negative 
arrivals  Fig.  3(c-d)  shows  that  Cases  3-4  have  smoother  and  more  consistent  arrival  tracking.  This 
difference  is  likely  due  to  the  varying  spatial  and  temporal  features  of  the  noise  field  generating  these 
coherent  wavefronts  (i.e.  individual  loud  ship  events  for  positive  arrivals  versus  more  diffuse  shipping 
noise  for  negative  arrivals).  Furthermore,  as  expected,  the  peak  SNR  is  highest  around  the  specific  day  or 

time  interval  that  is  used  to  construct  each  of  the  reference  covariance  matrices  ^  as  m^gd 

on  the  horizontal  axis  of  Fig.  3  with  a  red  line.  This  suggests  that  updating  the  reference  covariance 

matrices  ^ re-  ^  ^  during  the  span  of  the  total  observation  period  may  enhance  the  robustness  of  the 

STF  technique  for  noise-based  ocean  monitoring  (e.g.  by  using  a  different  reference  covariance  matrices 
for  each  day  of  the  experiment). 
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over  6  days  for  four  different  sets  of  steering  vectors  (a-d)  which  were  selected  from  the  same  Cases  1-4 
used  in  Fig.  2(f).  The  beamformer  output  was  computed  for  the  ensemble  averaged  cross-correlation 
waveforms  between  both  VLAs  obtained  using  a  short  moving  average  of  N=480  min  starting  at 
successive  minutes  L  over  6  days  (i.e.  1<L<8160).  The  index  L  thus  represents  the  beginning  of  the 
moving  time-window  along  the  “Slow  Time”  axis. 


The  total  cumulated  time  shift  over  6  days  was  found  to  be  positive  for  Cases  1-4,  and  the  shifts 
amounted  respectively  to  +1.1  ms,  +0.7  ms,  +1.1  ms,  and  +1.1  ms  (see  Fig.  4(a-d)).  Note  that  the 
cumulated  time  shift  has  the  same  sign  (positive)  when  tracking  coherent  wavefronts  occurring  at  positive 
time-delays  (Cases  1-2)  or  negative-time  delays  (Cases  3-4).  Thus,  this  time  shift  corresponds  to  an 
asymmetric  (or  non-reciprocal)  change  in  arrival  times  when  comparing  positive  and  negative  coherent 
wavefronts.  Flence,  this  overall  asymmetric  time-shift  is  likely  dominated  by  a  relative  clock  drift  of 
-200  ps/day  between  the  VLAs  autonomous  recording  systems  (which  causes  an  imperfect 
synchronization  of  the  noise  records)  as  its  trend  is  quasi-linear  over  six  days  and  most  likely  not  from  an 
actual  ocean  current  flowing  between  the  VLAs  [Sabra  et  al.,  2005b].  Furthermore,  other  environmental 
causes  (e.g.  variations  in  the  ocean  sound  speed  or  separation  distance  between  the  VLAs)  would  have 
otherwise  caused  symmetric  (or  reciprocal)  time-shifts  of  opposite  signs  with  respect  to  t=0  for  positive 
and  negative  coherent  wavefronts  [Sabra  et  al.,  2005b]. 

CONCLUSIONS 

The  results  of  this  research  effort  confirm  the  possibility  of  tracking  small  fluctuations  of  deterministic 
coherent  noise  arrivals  between  acoustic  arrays  in  a  challenging  coastal  environment  dominated  by  non¬ 
stationary  and  non-uniform  shipping  noise.  Further  studies  are  required  to  fully  assess  the  potential  of 
noise-based  ocean  monitoring  and  tomography  over  very  long  observation  periods  and  larger  array 
separation  distances. 

IMPACT 

It  is  conjectured  that  the  results  of  this  study  could  help  develop  a  totally  passive  means  for  monitoring 
the  ocean  environment  using  only  ambient  noise.  A  potential  scenario  benefiting  from  the  proposed 
methodology  might  include  long-term  deployment  of  ocean  sensing  systems  requiring  minimum  power 
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consumption,  covert  operations  in  hostile  settings,  or  coastal  deployments  where  active  sources  are 
limited  by  environmental  regulations. 
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